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Abstract 



The volume operator is an important kinematical quantity in the non-perturbative ap- 
proach to four-dimensional quantum gravity in the connection formulation. We give a general 
algorithm for computing its spectrum when acting on four-valent spin network states, evalu- 
ate some of the eigenvalue formulae explicitly, and discuss the role played by the Mandelstam 
constraints. 
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1 Introduction 



The volume operator has emerged as an important quantity in the kinematics of 3+1- 
dimensional quantum gravity in the loop representation. It is the quantum analogue of 
the classical volume function, measuring the volume of three-dimensional spatial regions. 
Although not an observable of the pure gravity theory (in the sense of commuting with 
the gravitational Hamiltonian), it has an immediate physical interpretation, and becomes a 
genuine observable if the spatial regions are defined intrinsically by additional matter fields, 
for example, the constant-value surfaces of a scalar field variable. Moreover, the volume is a 
relatively simple function of the dreibein variable E, of the canonically conjugate pair (A,E), 
where ^ G „4 is the (complex) S'C/(2)-Ashtekar connection on the three-space S. 

The significance of the volume operator for the quantum theory derives from the fact 
that its eigenfunctions span the kinematical quantum state space of functions on A/Q, the 
space of connections modulo gauge transformations, and that these functions are known, at 
least in principle [1]. They are the so-called spin network states, which can be expressed as 
appropriately (anti-)symmetrized linear combinations of certain Wilson loop functions (traces 
of holonomics of the connection variable A). This, as well as the spectral discreteness of the 
volume operator, was first pointed out in [2]. 

An obvious task at this stage is the actual computation of its spectrum. We will describe 
in the following a general algorithm for computing the volume spectrum on four-valent spin 
network states, that is, spin networks made up of Wilson loop states with no more than four 
line segments meeting at each loop intersection. The eigenvalues of several classes of such 
states will be computed explicitly. We have shown previously [3] that the spectrum in the 
three- valent case is identically zero, thereby correcting a computational error in [2], where 
a non-vanishing spectrum was derived. (Because of the algebraic structure of the volume 
operator it is clear a priori that loop states with valence less than three are annihilated.) 

Our calculations will take place in a lattice-regularized framework which we have been 
advocating elsewhere as an appropriate tool for approximating the quantum Hamiltonian 
dynamics of gravity in the loop approach [4]. However, to our understanding the form of the 
lattice operator coincides with that of the continuum formulations [5,6] (when restricted to 
the subset of states that can be realized on a three-dimensional cubic lattice), and therefore 
our results are equally valid in those cases. Since the discussion of the volume operator is 
largely insensitive to the signature of space-time and the particular form of the Hamiltonian, 
we will for simplicity work within the real, SU (2)-setting. We only remind the reader that 
there is a version of Ashtekar's gravity based on a real canonical variable pair {A,E), as 
recently discussed in [7] (where, however, the Hamiltonian assumes a more complicated form 
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than in the complex formulation). 

As part of the spectrum calculation, we will have to address the issue of over-completeness 
of the spin network states, which is a familiar feature of complete sets of gauge-invariant 
functions on a space A of gauge potentials. We will argue that it is necessary to eliminate 
this over-completeness in order to derive the correct spectrum for the volume operator. 

In the next section, we will recall the construction of the classical volume function and the 
discrctizcd quantum volume operator on the lattice. In Sec. 3, wc derive a general expression 
for the matrix elements of the volume operator on four-valent spin network states, which is 
then illustrated by explicitly calculating some of the spectral formulae in Sec.4. In Sec. 5 we 
summarize and discuss our results. 



2 Defining the volume operator 

In order to fix the notation, we will first summarize the main ingredients of the customary 
Hamiltonian lattice formulation for gauge theories [8] . Our lattice is a cubic NxNx A^-latticc, 
with periodic boundary conditions, i.e. the topology of a three-torus. The basic operators 
associated with each lattice link / are an SU (2)-link holonomy V, together with its inverse 
V~^, and a pair of canonical momentum operators pf and p~ , where i is an adjoint index. 
The operator pf{n,a) is based at the vertex n, and is associated with the link I oriented in 
the positive a-direction. By contrast, p~{n + la, a) is based at the vertex displaced by one 
lattice unit in the a-direction, and associated with the inverse link l~^{a) = l{—a). The wave 
functions are elements of ^iL^{SU{2),dg), with the product taken over all links, and dg is 
the Haar measure. The basic commutators are 



[VA''{n,a),Vc''{mM=0 

\pf{n,a),VA'-'{m,b)] = -| (5„,„(5-g TjA^Vs'^Cn, a), 

\p-{n,a),VA^{m,b)] = -'-5nmS^i,VA''{n,a)TiB^, (2-1) 

\pf{n,a),pf{m,b)] = ±iSnmS^i,eijkpt{n,a), 
\pf{n,a),pj{m,b)] = 0, 

where eijk are the structure constants of SU{2). In terms of an explicit parametrization by 
four complex parameters Oj, i = . . .3, a| = 1, the operators for a single link (n, a) are 
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given by 



^ \ —OLo + la-i an — la^ " /-^ * " 



= ^(ai^o - "o^i ± "3^2 =F "2^3), ^2 2) 

P2 = 2(^2^0 T 03^1 - 00^2 ± ai^a), 
^•3" = ^(«3<9o ± a25i =F ai32 - ao^s), 

where in the first line we have defined the three r-matrices. In the continuum theory, the 
classical expression for the volume of a spatial region 7?. C S is given by 



V(7^) = d^x Vd^ = ^ S-x yi^^^I^^i^f^^, (2.3) 

where are the dreibein variables introduced earlier (corresponding to the generalized 
electric fields in a gauge theoretic language) . Taking into account the continuum limit of the 
classical lattice variables, 



Vt^n^a)"-^ o'~E'l{n)^0{a\ (2.4) 
as the lattice spacing a goes to zero, we define the lattice analogue of (2.3) as 

Viatt = X] \/ ^'^^ (Pii'f^^^) +P7in, a)){p+{n,b) + pj {n,b)){p^{n,c) +p^{n,c))\. 

(2.5) 



For consistency, we have averaged over the momenta of both orientations"*^. The translation 
of this expression to the quantum theory is not well defined a priori, because of the presence 
of both the modulus and the square root. However, since both the p^ and therefore also the 
operators 



^ This choice does not alter the conclusions of [3], where the unaveraged operator was used. The eigenvalues 
reported in [3] are merely changed by a constant overall factor, e.g., C = ^ for four-valent intersections. 



Din) := ^eabce"^\pt{n,a)+Pi (n, a))(p+(n, 6) + p^. {n,b)){p+ {n,c) + (n,c)) (2.6) 



arc sclfadjoint, we may go to a Hilbert space basis of simultaneous eigenfunctions of all the 



through the square roots of the moduli of the eigenvalues of the D{n) in that basis. (Note that 
no operator ordering problem occurs in the definition of D{n).) As already mentioned in the 
introduction, the diagonalization of the volume operator is most easily achieved starting from 
a set of spin network states on the lattice. These are certain (anti-)symmetrized, real linear 
combinations of Wilson loops. (A Wilson loop on the lattice is a gauge-invariant function of 
the form TtV{li)V{l2)-.-V{lk), where j = h o I2 o ... o Ij^ is a closed loop of lattice links.) A 
spin network associates a positive "occupation number" with each lattice link, counting the 
number of (unoricnted) flux lines of basic spin-^ representations along the link, and also keeps 
track of the way in which those flux lines are contracted gauge-invariantly at the vertices (see 
[1] for more details). 

The operators D{n) have a particularly simple action on spin networks, because they do 
not change their support (in terms of the fliix line numbers). Thus only finite-dimensional 
rearrangements occur within each subset of states sharing the same occupation numbers, and 
the diagonalization of Viatt can be performed separately in these finite-dimensional eigcnspaces 
[2]. In this respect, the structure of the volume operator is much simpler than that of the 
Hamiltonian operator, which (at least on the lattice) changes the support of Wilson loops [4] . 
Since an operator D{n) acts only on links adjacent to the vertex n, and neighbouring Z)(n)'s 
commute, it is sufficient to study its action on spin networks locally around a single vertex. 
This will be the subject of the next section. 



3 Deriving the spectrum on four-valent spin networks 

Because of the cubic geometry of the lattice, the spin networks that can be defined on it 
are at most six-valent. 



D{n) and define the operator 




(2.7) 



n 
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Fig.l illustrates our labelling of the link directions meeting at a vertex n. One ingredient 
in describing a spin network locally around n is a 6-tuple j of integers ji > giving the 
occupation numbers (j'l, . . . , jg) of the links ((n, 1), (n, 2), (n, 3); {n, —1), (n, —2), (n, —3)) = 
((n, 1), (n, 2), (n, 3); (n — 1, 1), (n — 2, 2), (n — 3, 3)) intersecting at n. Since the flux lines are 
to be contracted at the vertex, their sum j := Yl^=i Ji an even integer. Next one has to 
specify how the j flux lines are joined pairwise at n to ensure gauge-invar iance. 

By convention we may join a flux line along the positive 1-direction, say, only to a flux 
line from one of the other five links, and not from the same link (i.e. we forbid "retracings" ) . 
This leads to a constraint on the occupation numbers: any ji has to be equal to or smaller 
than the sum of the remaining jk, i.e. ji < '^j^^ijk- Given a contraction of the flux lines 
at the vertex, the spin network consists of a weighted linear combination of Ji!j2'j3'i4b5'i6' 
Wilson loops corresponding to all possible permutations of flux lines associated with each 
of the six links. The weight factors are given by (— where P is the parity of the 
flux line permutation and N the number of closed loops in a multiple Wilson loop that is 
obtained by arbitrarily completing the local link configuration around the vertex n. The 
relative weights of the set of multiple Wilson loops thus obtained is independent of such an 
extension. To obtain a complete spin network state, this (anti-)symmetrization of course has 
to be performed around every vertex contained in the state. 

Since the diagonalization of the volume operator is algebraically rather involved, we will 
restrict ourselves to the simplest non-trivial case, namely that of spin networks with at most 
four-valent intersections. This is consistent since the volume operator maps the set of such 
states into itself. In a previous paper we have shown that spin networks are annihilated 

locally by D{n) at trivalent intersections, and more generally at intersections for which there 
exists only a single, unique contraction of flux lines [3]. From our calculations on four-valent 
spin networks below one recognizes this as part of a general pattern, namely that eigenvalues 
occur in pairs of opposite sign. 
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Without loss of generality we may restrict our attention to four-valent vertices with 
occupation numbers of the form j = {p, q,r; s, 0,0), with p>q>r>s>0. A moment's 
consideration reveals that the number of spin network states one may construct from this link 
configuration is {x + l)(a; + 2)/2, where x = Min{s, ^{—p + q + r + s)}. However, these spin 
networks are not all linearly independent, due to the existence of the so-called Mandelstam 
constraints (see, for example, [9]). It turns out that the number of linearly independent four- 
valent spin networks is x + 1. Thus it grows only linearly in x, whereas the total number of 
states is proportional to x'^. (It is a question of semantics whether by "spin networks" one 
means the full set of (anti-)symmetrized states as introduced above or only an independent, 
already orthogonalized basis set - we have been using it in the former sense.) 




An alternative way of parametrizing the occupation numbers is given by four integers 
{a, /?, 7, (5), where p = a + (3 + S,q = a + 'y,r = P + 'y, and s = 5 (see Fig.2; by drawing a 
etc. as connected pieces of incoming and outgoing flux lines we do not mean to indicate that 
the associated spin network states will share the same routings at the intersection; it merely 
gives us another way of labelling a flux line configuration.). Note that this transformation is 
invertible. The variable change is useful because of the following construction. Abandon for 
the moment the restrictions p > q > r > s, keeping however the inequalities p < {q + r + s), 
9 < ip + f + s), r < {p + q + s), s < + q + r) and the condition (p + q + r + s)mod 2 = 0. 
Given a fixed x > 1, we define the "fundamental link configuration" to be the one with 
minimal total occupation number (p + q + r + s). For this configuration we obviously have 
{p,q,r,s) = {x,x,x,x), or (a,/3, 7,5) = (0,0, x,x). Now observe that every allowed link 
configuration with the same x can be obtained by adding 4-tuplcs of non-negative integers 
(Aa, A/3, A7, A5) to (a, (3,^,6) = {0,0, x,x). The converse is not true, because adding such 
a quadruplet may change x. 
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Our aim is to derive general formulae for the eigenvalues of b{n) on independent spin 

networks for fixed x. This can be achieved in a series of steps. Starting from a link config- 
uration (a,/3, 7,(5) = (0,0, x,x), one extends it arbitrarily to obtain a multi-loop, and then 
determines the corresponding set of (x + + 2) /2 spin networks. A subset of x + 1 linearly 
independent states can be easily determined using, for example, the relations derived in [9]. 
Note that, due to the total (anti-)symmetry of the spin network states, loop configurations 
containing retracings on any of the links meeting at n can be set identically to zero. Next one 
computes the action of D{n) on the independent states, thus obtaining a {x + 1) x {x + 1)- 
matrix, which one then diagonalizes. It is equivalent and much easier to compute the action 
of D{n) on a single representative of the permutation equivalence class of Wilson loop states 
that make up a given spin network, and then check which permutation equivalence classes 
the resulting Wilson loop functions lie in. In all of this, one must remember to keep track of 
the weight factors (— A useful identity in evaluating the action of D{n) is 

e'^Hri)A^{rj)c''{Tk)E^ = 2 (S^S^S^ - S^SES^). (3.1) 

All of the above steps may be implemented using algebraic computing programs, like Mathe- 
matica. In the derivation of the eigenvalue formulae at constant x it is crucial to observe that 
it is sufficient to know the action of D{n) on a small number of link configurations (a, (3, 7, 5) 
(or rather their associated spin networks) "close" to the fundamental one, {0,0, x,x), as we 
will now proceed to explain. 

Let us adopt the shortcut described above for deriving the action of D{n) on spin networks 
by their action on Wilson loop representatives. Note that going from the link configuration 

{0,0, x,x) to {1,0, x,x), say, at the level of these Wilson loop representatives may be repre- 
sented by adding a single closed loop "containing a" (for an explicit example, see Sec. 4). It 
is then easy to repeat the steps outlined above, and determine the matrix elements of D{n). 
Call the matrices obtained in this way M(0,0,,t,,t) and A4{l,0,x,x) respectively. Using the 
explicit action of D{n) on Wilson loop states, one may then prove that any M{a,0,x,x), 
a = 2, 3, . . ., can be computed via 

M{a, 0, X, x) = M(0, 0,x,x)+a{ M(l, 0, x, x) - M(0, 0, x,x)). (3.2) 

Similar relations exist for other link configurations, and one can derive a general formula 
expressing the matrix representation of D{n) acting on the set of spin networks associated 
with an arbitrary link configuration (a, (5, 7, 5), 
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M^{a, A 7, (5) = (1 - 7 + x)(l - (5 + x) ( (1 - a)(l - (3) M,(0, 0, x, x) 
+ a{l - P)M^{l,0,x,x) + (1 - a)l3Mj;(0,l,x,x) + M^(l, 1, x, x) ) 
+ ( (7 - x)(l - (5 + + (1 - 7 + - x) ) ( (1 - a)(l - /3) M^(0, 0, a; + 1, x) 
+ a{l - (5) M^(l, 0, X + 1, x) + (1 - a)/3 M^(0, 1, x + 1, x) + a/? 1, x + 1, x) ) 



where we have introduced a subscript x to indicate that the matrix elements of M refer to 
sets of spin networks with fixed x. Formula (3.2) is of course a special case of (3.3). In the 
derivation of (3.3) we have taken into account that Ma;(0, 0, x + l,x) = Mr(0, 0, + 1), 
iV4(l,0,x + l,x) = Mr(l,0,x,x + 1) etc. 

Given this expression, it is straightforward to establish the explicit eigenvalue formulae 
for a given x. To illustrate our method, we will in the next section discuss the case x = 1 in 
some detail, and also give the eigenvalue formulae for x = 2 and 3. 



4 Computation of the spectrum for small x 

Recall first that, for x = 0, following [3], all eigenvalues of the volume operator vanish 
identically. Let us therefore turn to the case x = 1. The fundamental link configuration is 
given by (a,/3,7,(5) = (0,0,1,1) or, equivalently, {p,q,r,s) = (1,1,1,1), Fig.3a. 



r . . s 




p n Si S3 

Fig. 3a Fig. 3b Fig. 3c 

Now close up the external ends of the four links in some way, as shown schematically 
in Fig. 3b. There are then three spin networks, Si, S2 and S3, obtained by connecting the 
internal ends of the links meeting at the vertex n in various ways (Fig.3e). Since the links are 
singly occupied, no permutations have to be taken into account. The Mandelstam identity 
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for this set of spin networks is Si — S2 — Ss = 0. We choose as an independent set the states 
Si and 5'2. One then computes the matrix representation for D{n) on these states, which is 



given by Mi (0,0, 1,1) = ^ \^ ~1 J According to our reasoning of the previous section, 
we have to compute another seven matrices Mi. 

Consider next the hnk configuration (a,/3, 7, (5) = (1,0,1,1). Tliis means adding two 
links to the previous configuration, as ihustrated in Fig.4a and Fig.4b. The analogue of 
Fig.3c is shown in Fig.4c. 




Fig. 4a 




Fig.4b 





Wi 



W2 
Fig. 4c 




The three configurations Wi, W2 and W3 are now not spin networks themselves (because 
no (anti-)symmetrization has been performed), but can be considered as Wilson loop rep- 
resentatives of permutation equivalence classes, as explained earlier. For the corresponding 
spin networks S[, S2 and S'^, by construction S[ — S2 — S'^ = continues to hold, and one 



computes Mi (1,0, 1,1) = ^|^^ _2j^^ basis {S'i,S'2}. Proceeding similarly with the 
other relevant link configurations, one obtains for (3.3) 



Mi(a,/3,7,5) 



3i / (i + a)(l+/3) 



8 \il + a){P + ^ + 6) 



:i +/?)(« +7 +5) 

-(l + a)(l+/?) 



(4.1) 



from which one computes the eigenvalues of D{n) on link configurations (a,/?, 7, (5) as 



±0 V(l + 7 + ^)(1 + «)(1 + /?)(! + « + /9 + 7 + S)- 



(4.2) 



Note that we have used the identity (7 — 1){6 — 1) = 0, valid for x = 1. Going back to the 
notation (p, q, r, s) for link configurations, with p>q>r>s>0, one finds for the first few 
configurations the following eigenvalues: 
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(p,q,r,s) 


(1,1,1,1) 


(2,2,1,1) 


(3,2,2,1) 


(3,3,1,1) 


(3,3,3,1) 


(4,2,2,2) 


eigenvalues 


±|V3 




_l_3^/5 
^ 4 




_l_3^/3 
^ 2 


_l_3^/3 
^ 2 



Similarly, wc have computed the eigenvalues of D{n) for x = 2. For each link configura- 
tion there are six spin network states, out of which three are linearly independent. For the 
matrix M2(a, 7, 5), there is one zero-eigenvalue and a pair of non- vanishing ones. 



V-(2 + a)(2 + /3)(2 + a + /3) + (7 + 5 - 2)(a + /5 + 7 + (5)(2a/3 + 3a + 3/? + 4). 

4v2 

(4.3) 

In principle, formula (3.3) can be used also for the unreduced matrices, i.e. before the 
linearly dependent spin networks are eliminated. However, it is unfortunately not true that 
such states automatically have vanishing eigenvalues. Of course, if complex eigenvalues occur 
it is clear that they must correspond to spin networks that vanish modulo the Mandelstam 
constraints. But even non-zero, real eigenvalues occur for these spurious eigenvectors. (This 
happens, for example, in the case x = 3.) Again, we give a list of the first few eigenvalues 
corresponding to the matrix M2(a, /3, 7, (5): 



(p,q,r,s) 


(2,2,2,2) 


(3,3,2,2) 


(4,3,3,2) 


(4,4,2,2) 


(4,4,4,2) 


(5,3,3,3) 


eigenvalues 


^ 2 


_|_3VT3 
2V2 


2^/2 




±|V3 


±|V3 



For higher x, the algebra becomes progressively more complicated, but no problems 
occur in principle in determining the eigenvalues. Here are the eigenvalues of the matrix 
M^i^a, /3,j,6) for the case x = 3, which come in two pairs of opposite sign, namely, 



±1 sjh±^2, (4.4) 



where 
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h = 5fi{a + 2)(/3 + 2)u - 3^i{ii + 2) + 3{a + P + 6)(7 + 6) + 3(a/3 - 11) 

I2 = 16(/i + 2fiy'^{a + 1)(/^ + + 3)(/3 + 3) + 5(/Lt + 2)^1/2(50,2 + 5^2 ^ g^, ^ _ g^,^) 

+ 64(m + 2)^1/2 - 64(^ + 2y{a + 2)2(/3 + 2f 

+ I2{fi + 2)''^(10(a + 2)(/3 + 2) - 3(a + 2)^ - 3(/3 + 2)^)) 

+ 4(13Ai + 15)z/(a + 2)2(/3 + 2f - Qp'^{a + 2)(/3 + 2)(5a2 + 5/3^ + 48a + 48/3 + 244) 

+ 36zy2(2a^ + 17 + 44a + 2/3^ + 17/3^ + 44/3 + 72) - 94z/(a + 2f{f3 + 2f{a + (3 + 10) 

- 6z/(a + 2)(/3 + 2)(5a^ + 54q2 + 316a + 5/3^ + 54/3^ + 316/3 + 1152) 

+ 108//((a + 2)2(a + Af + {(3 + 2f{[3 + Af) + 36(a + Af{a + 2f + 36(/3 + 4)3(/3 + 2)^ 

+ 18(a + 4)(/3 + 4)(a + 2)2(/3 + 2f - 960(a + 2)(/3 + 2) 

+ 9(a + 2)(/3 + 2)(a + 4)(/3 + 4)(a2 - 26a + /J^ - 26/3 - 112) 

+ 6(a + 2)(/3 + 2)((a + 2)(-5a2 + 16a + 64) + (/3 + 2){-h0^ + 16/3 + 64)), 

(4.5) 

with the abbreviations /Lt = a + /3 + 7 + (5 and v = ^ + 5 — A. For small occupation numbers 
{p,q,r,s), some explicit eigenvalues are given by: 



(p^q^rjs) 


(3,3,3,3) 


(4,4,3,3) 


(5,4,4,3) 


(5,5,5,3) 


eigenvalues 


±|V19± 16 




±-^V33± V753 

4V2 


±^V33±27 



5 Summary and discussion 

We have derived a general formula for the representation matrix of the operator D{n) on 
sets of four-valent spin network states corresponding to a link configuration {p,q,r,s), with 
fixed X, where x = Min{s, ^{—p + q + r + s)}. This expression, formula (3.3), is given in 
terms of eight matrices whose entries have to be computed for each x. The spin network 
states are not independent, but obey certain linear relations, the Mandelstam constraints. 
The linearly dependent states have to be identified explicitly, otherwise one obtains spurious 
eigenvalues after diagonalization. At this moment, the only restriction on obtaining the full 
spectrum on four-valent spin networks is computing capacity. It may be possible that the 
formulae obtained above for x = 1, 2, 3 can be written in a more symmetric form that can be 
generalized to arbitrary x; so far we have not been able to do this. 
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We have observed that eigenvalues occur in pairs of opposite sign. Thus, for odd x + 1 
there is a zero-eigenvalue for an eigenstate that does not vanish modulo the Mandelstam con- 
straints. This means that there exist four-valent spin networks "without volume" (which by 
construction are non-planar; the planar ones arc all annihilated by V because of antisymme- 
try). Another consequence is that after taking the modulus, as is necessary for constructing 
the local volume operator V(n), all non-zero eigenvalues are (at least) two-fold degenerate. 

The generalization of our results to spin networks of higher valence is algebraically more 
complicated; in this case already the counting of spin networks is less straightforward. Still, 
there are no obvious obstructions to deriving analogues of our matrix formula (3.3). 

As for the geometric interpretation of our eigenvalue expressions, we have seen that their 
dependence on the occupation numbers {p, q, r, s) is not particularly simple, for example, 
they do not just depend on the total number of links, {p + q + r + s). One possibility is to 
try to extract from them a certain asymptotic behaviour, say, as the total number of links 
becomes very large. For example, for fixed x, one may look at spin networks corresponding 
to link configurations of the form {p,p,p,x), where p runs through all odd or even integers, 
depending on whether x is odd or even. For the cases studied in the last section, one finds 
the following asymptotic behaviour of the eigenvalues of D(n) for large p and to highest order 
in p: 



x=l: ~ '^~^P 

^ = 2: ~±^p2 (5.1) 
x = 3: ~±^(2±1)/ 

The "local volume", V(n), therefore increases linearly in p, and can become arbitrarily large 
for the classes of spin networks considered here. Such asymptotic relations may be useful in 
estimating the contributions of certain sectors of the Hilbert space of spin network states in 
numerical approximations. 

Acknowledgement. I am grateful to the members of the Max-Planck-Institut for their kind 
hospitality. 
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